Pairing Fluctuation Theory of Superconducting Properties in Underdoped to Overdoped Cuprates 
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We propose a theoretical description of the superconducting state of under- to overdoped cuprates, based on the 
short coherence length of these materials and the associated strong pairing fluctuations. The calculated Tc and 
the zero temperature excitation gap A(0), as a function of hole concentration x, are in semi-quantitative agree- 
ment with experiment. Although the ratio Tc / A{0) has a strong x dependence, different from the universal BCS 
value, and A(r) deviates significantly from the BCS prediction, we obtain, quite remarkably, quasi-universal 
behavior, for the normalized superfluid density ps(T) / ps{0) and the Josephson critical current Ic{T) / lc{0), 
as a function of T /Tc. While experiments on ps(T) are consistent with these results, future measurements on 
Ic{T) are needed to test this prediction. 



PACS numbers: 74.20.-z, 74.25.-q, 74.62.-c, 74.72.-h 

Pseudogap phenomena in the cuprates are of interest not 
only because of the associated unusual normal-state proper- 
ties, but more importantly because of the constraints which 
these phenomena impose on the nature of the superconductiv- 
ity and its associated high T^. Moreover, this superconduct- 
ing state presents an interesting challenge to theory: while 
the normal state is highly unconventional, the superconduct- 
ing phase exhibits some features of traditional BCS supercon- 
ductivity along with others which are strikingly different. 

Thus far, there is no consensus on a theory of cuprate su- 
perconductivity. Scenarios which address the pseudogap state 
below Tc can be distinguished by the character of the excita- 
tions responsible for destroying superconductivity. In the the- 
ory of Lee and Wen the destruction of the superconducting 
phase is associated with the excitation of the low-lying quasi- 
particles near the d-wave gap nodes. By contrast, Emery and 
Kivelson |Q] argue that the destruction of the superconductiv- 
ity is associated with low frequency, long wavelength phase 
fluctuations within a microscopically inhomogeneous model, 
based on one dimensional "stripes". 

In the present paper, we present an alternative scenario in 
which, along with the quasiparticles of the usual BCS theory, 
there are additionally incoherent (but not pre-formed) pair ex- 
citations of finite momentum q, which assist in the destruc- 
tion of superconductivity. This approach is based on a self- 
consistent treatment of the coupling of single particle and pair 
states. It represents a natural extension of BCS theory to the 
short coherence length (^) regime and provides a quantitative 
framework for addressing cuprate superconductivity. Here, 
we find a pronounced departure from BCS behavior in the un- 
derdoped limit which is continuously reduced with increasing 
hole concentration x. We derive a phase diagram for Tc and 
the zero temperature gap, A(0), as a function of x, which is 
in semi-quantitative agreement with (the anomalous) behav- 
ior observed in cuprate experiments, and we compute prop- 
erties of the associated superconducting state such as the su- 
perfluid density ps and Josephson critical cuiTent Ic- When 
these are plotted as ps{T) / ps{Q) and Ic{T)/ Ic{Q), as a func- 
tion of T/Tc, we deduce a quite remarkable, nearly universal 
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behavior for the entire range of x. 

As a simple model for the cuprate band structure, we con- 
sider a tight-binding, anisotropic dispersion ek = 2<||(2 — 
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cosfcj; — COS fey) + 2 tj_{l — cosfc^) — /i, 
is the hopping integral for the in-plane (out-of-plane) mo- 
tion and tj_ <C i|| [^. We assume that the fermions inter- 
act via an effective pairing interaction with d-wave symmetry 
Vk^k' = — IfflVkV'k' SO that (^k = ^ {cos kx — cos ky) . The 
present approach is built on previous work [Q-^] based on a 
particular diagrammatic theory, first introduced by Kadanoff 
and Martin and subsequently extended by Patton [||]. This 
approach can be used to describe the widely discussed BCS to 
Bose-Einstein cross-over problem [^, which has been associ- 
ated with small ^. The "pairing approximation" of Refs. [Q,^ 
leads to 



5= [l+5X(Q)]i(Q) 



(la) 



(lb) 



where ^{K) is the self-energy, x{Q) = 'Ek G{K) Go{Q - 
-^)'^k-q/2 P'^^^ susceptibility. Equations (1), along with 
the number equation n ~ 2 J^k G{K), self-consistently de- 
termine both the Green's function G{K) and the pair propa- 
gator, i.e., T-matrix t{Q). We use a four-vector notation, e.g., 
K EE (fe; z^), EE rE..,k and Go{K) ^ (to; - e^r\ 
We now show that these equations yield a natural extension of 
BCS theory to include incoherent pairs (labeled by pg), along 
with the usual quasiparticles and superconducting condensate 
(labeled by sc). 

We write the T-matrix and self-energy below Tc as t{Q) ~ 
tsc[Q) + tpgiQ), and J:{K) = ) + ^pg{K) . The con- 

densate contribution assumes the familiar BCS form tsc{Q) ~ 
—S{Q)Al^/T, where Age is the superconducting gap pa- 
rameter (and can be chosen to be real) and I]sc{K) = 
ft/ + Ek)- Inserting the above forms for the T-matrix 



into Eq. (|rb|), one obtains the gap equation 1 + .gx(0) = 0, 
as well as (for any non-zero Q), tpg{Q) = g/ (1 + gx{Q))- 
Note that because of the gap equation, tpg{Q) is highly 
peaked about the origin, with a divergence at Q = [|lO|]. 
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As a consequence, in evaluating the associated contribu- 
tion to the self-energy, the main contribution to the Q sum 
comes from this small Q divergent region so that T,pg{K) 
—Go{~K)Apg(pl^, where we have defined the pseudogap pa- 
rameter mJ 
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(2a) 



Thus, both Y.pg and the total self-energy S can be well approx- 
imated by a BCS-like form, i.e., ^{K) « tpl^/ {iujn + ek), 

where A = \J + ^pg magnitude of the total excita- 
tion gap, with the k dependence given by the d-wave function 
(ySk- Within the above approximations, the gap and number 
equations reduce to 
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(2b) 



(2c) 
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where the quasiparticle energy dispersion i?k = 
A^(p^)^/^ contains the full excitation gap A. 

The set of equations (^ can be used to determine the su- 
perconducting transition temperature Tc (where Asc = 0), 
and the temperature dependence of the various gap parame- 
ters. Eq. ( Pq ) contains the physics of the pair excitations, or 
pseudogap. The remaining two Eqs. ( ^b[ ) - ( ^ are analogous 
to their BCS counterparts but with a finite (as a result of non- 
zero Apg) excitation gap at Tc- 

It should be stressed that physical quantities which char- 
acterize the superconducting state depend on the pair and 
particle excitations, as well as condensate in different ways. 
The superfluid density can be written in terms of the London 
penetration depth as Ps,ab(T)/ps,ab{0) = [Aafc(0)/Aafc(T)]^, 
where 
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During the calculation special attention should be paid to lat- 
tice effects [ [l2[ ] and to the vertex correction (associated with 
the pseudogap self-energy) which enforces gauge invariance 
via the generalized Ward identity. This identity insures that 
Ps (X Ag^ and it vanishes identically at and above Tc- The 
prefactor A^^ = A^ — A^^ in Eq. (^ indicates that pairs (in 
addition to quasiparticles) serve to reduce the superfluid den- 
sity. 

In a related fashion, we address c-axis Josephson tunnel- 
ing between two identical high Tc superconductors. This 
situation is relevant to both break junction experiments [ pj| ] 
and to intrinsic Josephson tunneling Jl4| ] as well. An ex- 
pression for the Josephson critical current | p3| ] can be de- 
rived under the presumption that the tunnehng matrix element 
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FIG. 1. Phase diagram showing A(0) and as well as Apg{Tc) 
The inset shows experimental results from Ref. [p^]. 



iTkpP = iTol (5k||P|| + ll'il , where only the first (coherent) 
term contributes for a d-wave order parameter. 
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Equation (Q), like Eq. (g), differs from the usual BCS form 
(as well as that assumed by Lee and Wen [[l|,p^) in that the 
prefactor A^^ is no longer the total excitation gap A^. 

The remainder of this paper is directed towards understand- 
ing three experimental characteristics of the cuprates: (i) the 
phase diagram, (ii) the superfluid density and (iii) the Joseph- 
son critical current. 

(i) In order to generate physically realistic values of the var- 
ious energy scales, we make two assumptions: (1) We take 
g as doping-independent (which is not unreasonable in the 
absence of any more detailed information about the pairing 
mechanism) and (2) incorporate the effect of the Mott transi- 
tion at half filling, by introducing an x-dependence into the 
in-plane hopping matrix elements iy, as would be expected in 
the limit of strong on-site Coulomb interactions in a Hubbard 
model [|l7]| . Thus the hopping matrix element is renormalized 
as i|| (x) w to(l — n) = tox, where to is the matrix element 
in the absence of Coulomb effects. This x dependent energy 
scale is consistent with the requirement that the plasma fre- 
quency vanish at a; = 0. These assumptions leave us with one 
free parameter —g/Ato, for which we assign the value 0.15 to 
optimize the overall fit of the phase diagram to experiment. 
We take tj_/tn « 0.01 [|l|], and tQ sa 0.6 eV, which is reason- 



ably consistent with experimentally based estimates 1 19]. 

The results for Tc, obtained from Eqs. (^, as a function of 
X are plotted in Fig. [l|. Also indicated is the corresponding 
zero temperature excitation gap A(0) as well as the pseudo- 
gap Apg at Tc. These three quantities provide us, for use in 
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FIG. 2. Temperature dependence of the ab-plane inverse squared 
penetration depth. Main figure: from bottom to top are plotted for 
X = 0.25 (BCS limit, dot-dashed line), 0.2 (long-dashed), 0.155 
(dotted), 0.125 (dashed) and 0.05 (solid line). Lower inset: energy 
gaps as a function of T/Tc for x = 0.125. Upper inset: (A) the 
slope given by the low temperature expansion assuming a constant 



A^{0) 

(A-l-B) the sum of two contributions. 



Asc(T) = A(0); (B) the ratio 



at T/Tc = 0.2; and 



subsequent calculations, with energy scales which are in rea- 
sonable agreement with the data of Ref. |pj|], shown in the in- 
set. The temperature dependences of the energy gaps in Fig. 1 
are shown as the lower inset to Fig. 2, for a slightly under- 
doped case with x — 0.125. The relative size of Apg{Tc), 
compared to A(0), increases with decreasing x. In the highly 
overdoped limit this ratio approaches zero, and the BCS limit 
is recovered. This inset illustrates the general behavior as a 
function of T/Tc,: the excitation gap A is, generally, finite at 
Tc , the superconducting gap A^c is established at and below 
Tc, while the pseudogap Apg decreases to zero as T is reduced 
from Tc to 0. This last observation is consistent with general 



expectations for A 
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pg 
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It is important to stress, that our subsequent results for the 
superfluid density and Josephson current, need not be viewed 
as contingent on the detailed x-dependence used to derive the 
phase diagram. One can approach the calculations of these 
quantities by taking Tc{x) and the various gap parameters 
(shown in the inset) as phenomenological inputs, within the 
context of the present formalism. 

(ii) The superfluid density (normalized to its T = 
value), given by Eq. (^, is plotted in Fig. 2 as a function 
of T/Tc for several representative values of x, ranging from 
the highly over- to highly underdoped regimes. These plots 
clearly indicate a "quasi-universal" behavior with respect to x: 
Ps{T)/ ps{0) vs. T/Tc depends only slightly on x. Moreover, 
the shape of these curves follows closely that of the weak- 
coupling BCS theory. The, albeit, small variation with x is 
systematic, with the lowest value of x corresponding to the 
top curve. Recent experiments provide some preliminary ev- 
idence for this universal behavior IEOIeiIi. However, a firm 




FIG. 3. Temperature dependence of c-axis Josephson critical cur- 
rent for doping concentrations corresponding to the legends in Fig. ^ 
(main figure). The variation of A^;, and A as a function of T/Tc are 
plotted in the insets for the corresponding values of x. 



confirmation requires further experiments on a wider range of 
hole concentrations, from extreme under- to overdoped sam- 
ples [ p2|]. T his universal behavior appears surprising at first 
sight [|l[|l6|] because of the strong x dependence in the ra- 
tio Tc/A(0) (see Fig. 1). It should be noted that universal- 
ity would nof persist if ps{T) / ps{0) were plotted in terms of 
r/A(0). The nontrivial origin of this effect has a simple ex- 
planation within the present theory. At low to intermediate 
temperatures, it follows from Eq. (§) that A^^(0)/A^f,(T) = 
(AL(r)/A2(0)) (1 - AT/Tc + 0[iT/Tc)']) ^1-[A + 
B(T)] {T/Tc), where terms of order {T/Tcf and higher 
have been neglected. Here A ?,2\/2\n2 {e^ / c^) t\\ 
X (Tc/A(0)) represents the standard contribution to the lin- 
ear T dependence of Ps{T). The new term B{T) = 
(Tc/T) Alg (T) / A^ (0) derives from the pseudogap contribu- 
tion and has a weaker than linear T dependence (as can be 
inferred from the lower inset in Fig. 3) [p3[|. For the purposes 
of illustration, these two terms are plotted in the upper inset 
of Fig. 2 at T/Tc = 0.2. Note that the effective (negative) 
"slope" A + B is relatively x independent over the physical 
range of hole concentrations. Physically, the terms A and B 
are associated with two compensating contributions, arising 
from the quasiparticle and pair excitations, respectively, so 
that quasi-universal behavior results at low T. It can be shown 
that the same compensating effect obtains all the way to Tc, 
as is exhibited in Fig. 2. Thus, the destruction of the super- 
conducting state comes predominantly from pair excitations 
at low X, and quasiparticle excitations at high x. 

(iii) Finally, as plotted in Fig. 3, we obtain from Eq. (Q), 
similarly, unexpected quasi-universal behavior for the normal- 
ized c-axis Josephson critical current for the same wide range 
of X as in Fig. 2. This behavior is in contrast to the strongly 
X dependent quasiparticle tunneling characteristics which can 
be inferred from the temperature dependent excitation gap 
plotted in the upper inset of Fig. 3. The origin of this univer- 
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sality is essentially the same as that for pa, deriving from two 
compensating contributions. At this time, there do not appear 
to be detailed studies of Ic{T) as a function of x, although fu- 
ture measurements will, ultimately, be able to determine this 
quantity. In these future experiments the quasiparticle tunnel- 
ing characteristics should be simultaneously measured, along 
with Ic{T), so that direct comparison can be made to the exci- 
tation gap; in this way, the predictions indicated in Fig. 3 and 
its upper inset can be tested. Indications, thus far [ pj| , p4| ], are 
that this tunneling excitation gap coincides rather well with 
values obtained from photoemission data (see Fig. 1). 

In summary, in this paper we have proposed a scenario for 
the superconducting state of the cuprates. This state evolves 
continuously with hole doping x, exhibiting unusual features 
at low X (associated with a large excitation gap at Tc) and 
manifesting the more conventional features of BCS theory at 
high X. In this scenario the pseudogap state is associated with 
pair excitations, which act in concert with the usual quasipar- 
ticles. Despite the fact that the underdoped cuprates exhibit 
features inconsistent with BCS theory (rc/A(0) is strongly x 
dependent and A is finite at and above Tc) we deduce an in- 
teresting quasi-universality of the normalized ps and Ic as a 
function of T/Tc- In these plots the over- and under-doped 
systems essentially appear indistinguishable. Current experi- 
ments lend support to this universality in ps, although a wider 
range of hole concentrations will need to be addressed, along 
with future systematic measurements of Ic- 
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